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Photosynthetic complexes improve the transfer of excitation energy from peripheral antennas to 
reaction centers in several ways. In particular, a downward energy funnel can direct excitons in 
the right direction, while coherent excitonic delocalization can enhance transfer rates through the 
cooperative phenomenon of supertransfer. However, isolating the role of purely coherent effects 
is difficult because any change to the delocalization also changes the energy landscape. Here, we 
show that the relative importance of the two processes can be determined by comparing the natural 
light-harvesting apparatus with counterfactual models in which the delocalization and the energy 
landscape are altered. Applied to the example of purple bacteria, our approach shows that although 
supertransfer does enhance the rates somewhat, the energetic funnelling plays the decisive role. 
Because delocalization has a minor role (and is sometimes detrimental), it is most likely not adaptive, 
being a side-effect of the dense chlorophyll packing that evolved to increase light absorption per 
reaction center. 


Photosynthetic organisms harvest light using antenna 
complexes containing many chlorophyll molecules [I] . The 
energy collected by the antennas is then transmitted, 
through excitonic energy transfer (EET) [2], to a reaction 
center (RC), where it drives the first chemical reactions 
of photosynthesis. The thorough study of EET in photo¬ 
synthetic antennas has been motivated, in part, by the 
prospect of learning how to design more efficient artificial 
light-harvesting devices mill. 

It has long been recognized that excitons in many pho¬ 
tosynthetic complexes are directed toward the RC ener¬ 
getically: if the antennas lie higher in energy than the 
RC, the excitons can spontaneously funnel to the RC. A 
more recent discovery is that coherent mechanisms can 
also enhance light-harvesting efficiency. In particular, ex¬ 
citonic eigenstates may be localized or delocalized over a 
number of molecules, depending on the strength of their 
couplings laiMH]. Delocalization—i.e., coherence in the 
site basis—makes the aggregate behave differently than 
a single chlorophyll molecule, and phenomena such as 
superradiance [H [10] , superabsorption m, and super¬ 
transfer HHE] can occur in densely packed aggregates. 
Specifically, supertransfer occurs when delocalization in 
the donor and/or the acceptor enhances the rate of the 
(incoherent) EET between them. 

The presence of both funnelling and supertransfer sug¬ 
gests that their contributions to the efficiency could be 
quantified. However, the two effects are too closely related 
for such a separation to be easily carried out; in partic¬ 
ular, a change to the extent of delocalization requires 
changing excitonic couplings, which also determine the 
energy landscape. In other words, because delocalization 
and the energy landscape are intimately connected, it is 
not sufficient to alter one property to see what happens 
to the efficiency, because doing so also alters the other 
property as well. 


Here, we show that the roles of the two processes can be 
separated by constructing counterfactual light-harvesting 
complexes that reflect plausible evolutionary alternatives. 
In particular, screening thousands of complexes with vary¬ 
ing energetic landscapes and extents of delocalization al¬ 
lows the effects of changing one property to be examined 
while keeping the other as constant as possible. 

We describe this approach through its application to 
the strikingly symmetric antenna complexes of purple 
bacteria [16], which feature tightly packed bacteriochloro- 
phylls and considerable excitonic delocalization [nHiH]. 
This delocalization is known to give rise to supertransfer, 
in particular for EET within the LH2 complex [29H3T] . 
In principle, supertransfer between different complexes 
should also occur, and although it has not been studied as 
well as intra-complex supertransfer, it has been proposed 
to explain the high light-harvesting efficiency [IlEz]. 

We show that although supertransfer is present in pur¬ 
ple bacteria, it is not essential for efficient light harvesting. 
Eor example, EET efficiency can remain roughly as high 
even if half of the pigments are removed, weakening the 
couplings and localizing the states. When delocalization 
is removed from only some of the complexes, the efficiency 
can change drastically, but these changes are almost en¬ 
tirely due to shifts in energy levels; even in the worst cases, 
modest modifications of the site energies can restore the 
high efficiency. Indeed, in the presence of a strong fun¬ 
nel, the efficiency can be high regardless of delocalization. 
This suggests that the evolutionary advantage of densely 
packed chlorophylls is that they enhance the absorption 
cross-section per RC, while the delocalization is merely a 
side-effect of the dense packing. In evolutionary language, 
delocalization is a spandrel, not an adaptation [32] . 
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MODEL 


Structure 


We study the photosynthetic apparatus of the purple 
bacterium Rhodohacter sphaeroides^ which consists of re¬ 
action centers surrounded by two types of transmembrane 
antenna complexes, LHl and LH2 m- Although all 
complexes can absorb light directly, the dominant energy- 
transfer pathway is LH2 ^ LHl ^ RC. The complexes 
can arrange themselves in many ways [33l [34] ; as a rep¬ 
resentative example we consider an array of ten evenly 
spaced LH2s surrounding the core LHl-RC complex, as 
shown in Fig. 

The coordinates of the pigments in Fig. ^ are taken 
from the available crystal structures [35l|36]. Each LHl 
complex consists of 56 bacteriochlorophyll a molecules 
(BChl) with an average nearest-neighbor Mg-Mg distance 
of 8.5 A and immediately encircling two RCs [36]. These 
BChls absorb at 875 nm and are thus referred to as the 
B875 aggregate. Each RC comprizes four BChls, two in 
the tightly-coupled special pair and two accessory ones. 
The absorption spectrum of the RC has peaks at 865 nm 
and 802 nm, corresponding to the special pair and the 
accessory BChls, respectively. 

LH2 complexes surrounding the LHl-RC core include 
two rings each, one with 9 BChls (the B800 subunit) 
and the other 18 (B850) [35]. We only consider the 
B850 ring because B800 excitons are transferred to B850 
efficiently and quickly (700 fs [37]) and efficiently using 
supertransfer [291I3T] . making their ultimate fate nearly 
identical to that of excitons starting at B850. The average 
Mg-Mg distance between BChls in B850 is 9.0 A. 

The distances between complexes have been determined 
by atomic force microscopy (AFM) [33l|34], and our model 
uses the most common center-to-center spacing between 
LH2s, 75 A. This gives a distance of 22.5 A between 
nearest Mg atoms in different LH2s, which we also used 
to set the separation between the LH2s and LHl. 

We refer to the geometry just described as the nat¬ 
ural geometry (N), to distinguish it from the trimmed 
geometry (T) (Eig. Ht), in which every other BChl was 
removed to reduce nearest-neighbor couplings and encour¬ 
age exciton localization (see below). Eor example, the 
distance between BChls in B850 roughly doubles upon 
trimming, becoming comparable to that in the B800 ring, 
where delocalization is known to be minor. To trim each 
RC, out of the four BChls, we kept one of the special 
pair and the less-strongly coupled accessory BChl. We 
also considered cases where only some of the complexes 
(LH2, LHl, and/or RC) were trimmed. The extent of 
the trimming (keeping every second BChl) is not critical 
to our argument: Supplementary Table 4 shows similar 
results if only every third BChl is kept. 
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Figure 1. Model of the photosynthetic apparatus of 
Rh. sphaeroides, including the reaction center surrounded by 
the antenna complexes LHl and LH2 (only B850 subunit 
shown). Drawn using VMD [38] . 


Excitonic couplings 

Strong coupling between neighboring BChls in the nat¬ 
ural geometry leads to exciton delocalization [2 [IHS]. 
Within each complex (LH2, LHl, or RC), the excitonic 
states are the eigenstates of a Frenkel Hamiltonian [2 
that, in the weak-light, single-exciton regime, reads 

if = ^i;,|*)(i| + ^U,{|z)0-| + b-)(z|), (1) 

i i<j 

where \i) are the site states, Ei are the site energies, and 
Vij are the couplings. 

The site energies and couplings have been calculated by 
numerous groups [niiTiiisiEiiio], with the final values 
varying widely due to different electronic-structure meth¬ 
ods and assumptions about the molecular environment. 
For instance, values between 238 cm“^ and 806 cm“^ have 
been used for the coupling between the la and ip BChls 
in LH2 [40]. To confirm that our conclusions are general 
and not sensitive to the details of the models, we carried 
out the entire study using two sets of simulation parame¬ 
ters, which we call S and R because they are inspired by 
the approaches taken in Schulten’s and Renger’s groups 
(Table . Some differences between the two approaches 
are worth noting. 

Schulten obtains nearest-neighbor couplings from 
electronic-structure theory (Table and more distant 
ones using the point-dipole approximation (PDA) [TTIIT^ 
\2T\ [244[27] . Although the PDA is accurate at large 
separations, it can fail even at moderate separations 
of tens of angstroms, enough to encompass many non- 
nearest-neighbors [414(44] . By contrast, Renger’s TrEsp 
method [45] reproduces quantum-chemical couplings with 
an accuracy better than the unavoidable errors in almost 
all cases [44], which is why we use it for all couplings 
in R, except for the special pair, which is too close for 
TrEsp [46] . 
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Parameter set 





S' 

R 


Site 

LH2 

Eo, 

12 458 [21] 

12 m 

a 

energies 


A 

12 654 [m 

12 274 

! 

(cm“^) 

LHl 

Flhi 

12121 [E] 

11701 



RC 

Ep, 

12 180 [26] 

1199^ 

ir 



Fp2 

12 080 [21] 

1199E 

ir 



Ebi 

12 500 [21] 

12 47c 

ir 



Eb2 

12 530 [21] 

12 47c 

?F 

Nearest- 

LH2 

Vlaip 

363 [21] 

_1 

D 


neighbor 


Vip2a 

320 [21] 

_i 

7 


couplings 

LHl 

alternating 

300, 233 [21] 

_ i 

7 


(cm“^) 

RC 

RpiP2 

500 [21] 

418 [E] 



FPlBi, FP 2 B 2 

-50 [21] 

_] 

D 




Vp^B2 ? LP 2 B 1 

-60 [21] 

_i 
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More 


Method 

point dipole 

TrE^ 

FI 

distant 


Sr 

f] .rn 

1.2 



couplings 


m(D) 

rT.7^ 

6.1 




^ Chosen to align brightest state with absorption maximum. 

^ Computed in the same manner as the distant couplings. 

In S', we follow Ref. |21| in using alternating couplings even for 

the non-circular model of LHl. 

Using the transition charges in m 
® Corresponding to C = 348 000 cm ^ in | 47| . 

^ Corresponding to f = 0.8 in m 

Table 1. Site energies and couplings. Abbreviations: ol and 
(3: alternating BChls in the B850 unit of LH2; Pi, P 2 : special 
pair BChls in RC; Bi, B 2 : accessory BChls in RC. 

The largest differences between S and R are in the 
relative permittivity Sr and the transition dipole moment 
/i of each BChl, both of which enter coupling calcula¬ 
tions (whether PDA or TrEsp) through the ratio C = 
/i^/dTT^o^r- Schulten sets C so that the PDA reproduces 
the quantum-chemical coupling (in vacuum) between a 
particular pair of BChls [1711171119]. This can lead to er¬ 
rors because the long-range couplings now depend on the 
quantum-chemical method used for the short-range cou¬ 
pling between the calibration pair, as well as on the choice 
of that pair. Consequently, published values of C have in¬ 
cluded 116 000 [49], 146 798 [26|, 170 342 [20], 348 000 [47], 
and 519 310 cm“^ [17]. We use = 2 and /i = 11.75 D, 
corresponding to C = 348 000 cm“ 1 EniiiiiiT], while 
noting the unrealistically large dipole moment. 

Renger has argued that coupling calculations should 
use the measured value of fi miEn], estimated at 6.1 D in 
vacuum and 7.3 D in a medium with = 2 (the effective 
permittivity of protein complexes m) [52]. Calculations 
with molecules in cavities surrounded by the dielectric 
lead to an effective Sr from 1.25 to 1.67 m- The rec¬ 
ommended combination of fi = 6.1 D and Sr = 1.25 [48] 
corresponds to C = 149 000 A^ cm“^. (Alternatively, for¬ 
going the cavity model and using Sr = 2 and /i = 7.3 D 
gives a similar value, C = 133 000 A^ cm“^.) The large 
discrepancy between S and R on the value of C has a 


substantial influence on inter-complex FRET rates, which 
are proportional to . 

For the S parameter set, site energies are taken from 
Schulten’s papers (see Table [^, and R energies are chosen 
so that the brightest states of each complex is aligned 
with the absorption maximum. 

We neglect disorder in site energies, because its effect 
on the inter-complex transfer rates is small (about 20% for 
LH2^LH2 [24]) compared to that caused by differences 
between S and R parameters. The resulting difference in 
the efficiency would be even less, as is the case for the 
cyanobacterial photosystem I, where disorder changes the 
efficiency by about 1% m- 


Energy transfer rates 

The two geometries—natural and trimmed—differ in 
the inter-pigment couplings, which affects the excitonic 
states. In the T geometry, we assume that the coupling 
between the pigments is much weaker than the coupling 
of the pigments and their environment. In that case, ex- 
citons can be thought of as localized on individual sites, 
and EET is described using Forster’s theory of resonant 
energy transfer (FRET) [2]. By contrast, in the N ge¬ 
ometry, we assume that the inter-pigment coupling is 
much stronger than the pigment-bath coupling, leading to 
delocalized excitons being the better theoretical descrip¬ 
tion. For every aggregate (i.e., each LH2, LHl, and RC), 
we construct the excitonic states —the eigenstates 

of that aggregate’s Hamiltonian (Eq. —with dynamics 
described using Redfield theory [2]. 

The two approaches—full localization and full 
delocalization—are two ends of a spectrum. The reality— 
partially localized eigenstates, depending on the details 
of the system-bath interaction—is somewhere in between. 
Nevertheless, we focus on these two extremes because 
they serve to isolate the effect of delocalization on effi¬ 
ciency, which would persist in a diminished form even if 
the delocalization were only partial. 

In all cases, the coupling between different aggregates 
is weak and is described by FRET. However, delocaliza¬ 
tion within the aggregates can affect the inter-aggregate 
FRET rate, potentially leading to supertransfer [HHH], 
an effect proposed to be important for the efficiency of 
purple-bacterial light harvesting dlllT]. For example, 
if the donor contains two pigments and the exciton is 
localized on either with probability 1/2, each pigment 
contributes equally to the total EET rate to the accep¬ 
tor, k\oc oc + ^|/^D 2 /^Ap. But if the exciton is 

delocalized over the donor sites, the effective transition 
dipole of the donor is a linear combination of molecular 
transition dipoles. In the best case, with the exciton in 
the bright state {\Di) |D 2 ))/\/ 2 , the transfer rate is 

doubled, fcdeloc OC = 2fcioc if MDi = MDa- 
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Figure 2. Details of the model, (a) In the natural geometry (N), all the pigments are included. Due to strong intra-complex 
couplings, excitations are delocalized and inter-complex energy transfer is described by generalized FRET, which transfers 
population from donor eigenstates to acceptor eigenstates, (b) In the trimmed geometry (T), every second pigment is removed 
to weaken the couplings and localize the excitons, so that energy transfer is described by ordinary, site-to-site FRET. Only the 
larger FRET rates are shown, (c) Excitonic energy levels in the natural geometry {S parameter set). Brighter states are drawn 
thicker. Examples of each process included in our model are shown: energy transfer optical pumping radiative 

relaxation non-radiative relaxation internal conversion (k^^), and charge separation in the RC (k^^). (d) In the 

trimmed geometry, the site energies are not split by the excitonic couplings. The energy levels in the R parameter set are shown 
in Supplementary Eig. I. 


The approximation that the overall transition dipole 
of the donor eigenstate couples to the acceptor is only 
valid at large separations. Otherwise, the acceptor sees 
individual donor sites and not simply a super-molecule, 
a situation correctly described using generalized FRET 
(gFRET) [T3l |29l [53l |54] , which gives the coupling 

between a donor state 'ijj and an acceptor state (j) as 

= ( 2 ) 

hi 


where cf (c^) are the site-basis coefficients of the exciton 
states of the donor (acceptor) and Vij is the coupling 
between site i of the donor and site j of the acceptor. 
To distinguish ordinary, site-to-site FRET used in the T 
geometry from the gFRET in the N geometry, we refer 
to the former as oFRET. 

Once the coupling is known, the FRET rate is [2] 


k 


ET 

nm 


27r 

h 


Jn 


(3) 


where the indices m (donor) and n (acceptor) refer to 
either sites i and j in oFRET or eigenstates and (j) in 


gFRET. Jnm = / Lm{E)In{E) dE is the overlap between 
the normalized fluorescence spectrum Lm of m and the 
normalized absorption spectrum In of n. 

The spectra Lm and In can be calculated using mul- 
tichromophoric FRET theory, given a detailed model of 
each pigment’s environment (SU [55l [56] . For example, a 
gFRET calculation of the LH2 —> LH2 transfer rate us¬ 
ing an Ohmic spectral density [56] agrees with the same 
calculation with the exact HEOM method m- However, 
transfer rates are sensitive to changes in the spectral den¬ 
sity, meaning that substantial errors are introduced by 
the simplistic assumption of an Ohmic spectral density. 
Because we doubt that currently available spectral densi¬ 
ties are accurate enough to justify computing Lm and 
we parametrize Jnm using experimental data. 

We make the simplest possible choice, taking both 
Lm{E) and Im{E) to be normalized Gaussian func¬ 
tions centered at with the standard deviation a = 
250 cm“^ chosen to reproduce the width of the Qy band 
of BChl a in solution [Tj. With Gaussian spectra, Jnm 
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becomes 


Jnm = e ®"™/^'^"/\/47r(T2, (4) 

where Enm is the energy difference. Neglecting the Stokes 
shift between and Im would imply Jnm = Jmm caus¬ 
ing Eq. [^to break detailed balance. We correct this by 
using Eq. only for energetically downward transitions 
{Em > En), otherwise taking , 

where Tb = 300 K is the ambient temperature. 

Transfer rates obtained using this approach agree with 
those obtained by Schulten’s group. Eor LH2^LH2 
transfer, with the center-to-center distance increased to 
85 A, we obtain a transfer time of 13 ps, close to the 9.5 ps 
obtained by Schulten’s group using gERET with their 
chosen Ohmic spectral density m- 


Relaxation and optical pumping 

A complete model of EET also includes exciton loss 
(recombination) and creation by optical pumping. 

Excitons can recombine radiatively or non-radiatively. 
We incorporate radiative relaxation of site/exciton m 
by adding the rate of spontaneous emission = 

kf^\limg/l^()?{Erng/E()Y, where k^^ = (16.6ns)“^ is the 
radiative decay rate of BChl in solution [7], 
is the ratio of the oscillator strength of the transition to 
that of the single BChl, and (Emg/Eo)^ is the (small) 
correction to spontaneous emission as the energy of the 
transition changes with respect to the transition energy 
of BChl in solution, Eq = hc/(770nm). We include a 
non-radiative recombination rate of k^^ = (lns)“^ for 
all m [21]. 

In the natural geometry there is also intra-aggregate re¬ 
laxation (internal conversion) among the excitonic states, 
occurring on a sub-picosecond timescale m- We incor¬ 
porate it by assuming, within each aggregate (LH2, LHl, 
or RC), internal conversion from higher-energy states to 
lower ones at a rate k]^ = (100fs)“^, while the energeti¬ 
cally uphill rates are included by detailed balance. The 
speed of internal conversion relative to other processes 
means that each aggregate will be close to a Boltzmann 
state, as was confirmed for LH2 [24]. 

In natural light, the incoherent light populates system 
eigenstates and not their superpositions [58H6Q] . giving 
rise to a steady state In principle, these are vi- 

bronic eigenstates of the entire light-harvesting apparatus, 
which raises two considerations. First, eigenstates of the 
entire apparatus will include superpositions of different 
complexes. However, inter-complex delocalization is often 
neglected because it is destroyed by dynamic localization 
faster than other relevant timescales. This is especially 
true when the complexes are not energetically resonant 
and when the system-bath coupling is stronger than the 
inter-complex coupling, both conditions that apply here. 


Furthermore, although inter-complex delocalization may 
increase the absorption of the RC somewhat [63], that 
would modify the efficiency only slightly because most 
light is absorbed by the antenna complexes and not the 
RC. Second, the steady-state density matrix of each com¬ 
plex may not be diagonal in the electronic basis due to 
the system-bath coupling [niiMiiM]- However, our intra¬ 
complex Redfield treatment assumes weak system-bath 
coupling, meaning that off-diagonal elements will be small. 
We therefore neglect this correction because its influence 
is likely to be smaller than, say, the difference between 
S and R parameters. Therefore, we assume that inco¬ 
herent light populates the sites in the T geometry and 
the excitonic states of individual aggregates in the N ge¬ 
ometry. For site/exciton m, the optical pumping rate is 
k^ = k^{^n{Em)^ where n{Em) = — 1)“^ is 

the mean photon number at that energy and Tr = 5780 K 
is the effective black-body temperature of solar radiation. 

The final ingredient is the assumption that excitons 
in the RC can drive charge separation, at a rate k!^^ = 
k^^ = (3ps)“^, m G RC [HEl]. 


Master equation and efficiency 

With all the rates described above, and because there 
are no coherences in incoherent sunlight, the dynamics 
of the system can be described using a Pauli master 
equation. 


P = 


(5) 


where p is the vector containing the populations Pm of all 
the sites (T geometry) or excitonic states (N geometry), 
along with the population pg of the ground state, while 
the rate matrix K includes all the rates listed above (and 
summarized in Fig. [^-d). 


_ ^et 


uRR 


UNR 


ac 


k: 


OP 


k, 


cs 


(for n ^ m)^ 




— - ^ K 

n^m 


(6) 

(7) 


We define the efficiency as the quantum yield of charge 
separation, i.e., the probability that a photon absorbed by 
any of the complexes eventually drives charge separation 
in the RC. Because energy can be lost along the way, this 
is not a thermodynamic efficiency. 

Because incoherent excitation is stationary, the molecu¬ 
lar ensemble will be at steady state, [TSJ [61] [62] . Since 
pSS _ ^pSS _ pSS easily found as the unique 

eigenvector of K with eigenvalue zero. The efficiency is 
then the rate of charge separation in the RC divided by 
the rate at which excitons are created. 


T] = 


uCS y- ^SS 

^ Z^mGRC Pm 

„SS V A’OP ’ 

J^g Z-^m ^mg 


( 8 ) 
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■ Original energies, S 

■ Original energies, R 
□ Optimized energies, S 

■ Optimized energies, R 


Geometry 


Figure 3. Energy transfer efficiency for various geometries. The three letters indicate whether the natural (N) or trimmed (T) 
geometry was used, respectively, for LH2, LHl, and RC. The efficiencies were calculated using the energies in Tableas well as 
with energies optimized for maximum efficiency, and using both parameter sets S and R. The high efficiency in many cases with 
trimmed complexes indicates that delocalization is not essential for high efficiency. The large improvements when the energies 
are optimized indicates that energy-level alignment is more important for the efficiency than the effects due to delocalization. 


where is the steady-state population of the ground 
state. The in the denominator ensures rj is intensity- 
independent as long as the light is weak and ^ 1 . 

Eq. can be compared to the efficiency used by Schul- 
ten’s group, 77 = where is a uni¬ 

form distribution over the RC sites and is the initial 
distribution EOlllI]. This definition is not conceptually 
suited to steady-state light harvesting, where there is no 
“initial” state, although the two approaches are equivalent 
if a suitable “initial” state is chosen [65]. Because the 
initial state can influence the efficiency substantially m, 
it is important to model the light absorption and not 
assume an initially localized state mi HU or a uniform 
distribution over many [49] . 

RESULTS AND DISCUSSION 
Original model 

The overall efficiencies of the various cases are sum¬ 
marized in Fig. with details, including inter-complex 
transfer rates, in Supplementary Table 1. The results ob¬ 
tained using the S and R parameters can be quite different 
due to the weaker long-range coupling in R, confirming 
our concern about the difficulty of precise calculations 
and the need to focus on general trends. 

Figure shows that delocalization is not necessary for 
high efficiency. In particular, in the S parameters (but not 
R —see below) the trimmed geometry TTT is as efficient 
as the natural geometry NNN. In addition, delocalization 
in LH 2 always diminishes the efficiency: with both S 
and R parameters, tjtxy > ^nxy for ah X, Y G {N, T}, 
disproving the hypothesis that delocalization in LH2 ben¬ 
efits the efficiency through supertransfer. By contrast, 
delocalization in the RC is always beneficial, while in LHl 


delocalization benefits the efficiency if and only if there 
is delocalization in the RC. 

Whether delocalization is beneficial or deleterious can 
be explained using Eq. Delocalization affects the EET 
rate by altering both (potentially yielding super¬ 

transfer) and the excitonic energy gaps (through the Davy¬ 
dov splittings). Indeed, the spectral overlap is exponen¬ 
tially sensitive to excitonic splittings (Eq.[^, unlike super¬ 
transfer, which is at best an improvement by a constant 
factor equal to the number of sites. 

Delocalization in LH2 is deleterious because it weak¬ 
ens the energy funnel into LHl. LH2 sites are higher in 
energy than LHl sites, ensuring that backward transfer 
LHl ^ LH2 is suppressed. If the excitons are delocalized, 
the bright state of LH2 is lowered in energy by the Davy¬ 
dov splitting, bringing it closer to LHl. Although this 
encourages forward transfer, it increases backward trans¬ 
fer even more because of detailed balance. Furthermore, 
because there are 10 LH2s and only one LHl, there are 10 
times as many pathways back to LH 2 . Consequently, the 
highest efficiencies occur when the forward rate is much 
greater than the backward rate (see Supplementary Table 
1). In particular, the main pathway between LH2 and 
LHl is through their brightest eigenstates, which have a 
gap of about 350 cm“^. If LH2 is trimmed, EET occurs 
through the sites, whose energy is about 1000 cm“^ above 
the brightest state of LHl. The increased gap makes back¬ 
ward EET exceedingly slow, while keeping the forward 
rate adequate for high efficiency. 

The analysis in the RC is analogous. Its site energies 
are higher than in LHl, which would inhibit efficient 
transport. Including the excitonic couplings splits the 
levels, brining the lower one closer to LHl in energy and 
increasing the rate of forward EET [27]. The efficiency 
is worst when LHl is delocalized and RC is not (NNT 
and TNT), because that maximizes the energy difference 
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^LH 1 -^RC (1000 cm Elhi-Hrc (1000 cm Elhi-Hrc (1000 cm 

Figure 4. Efficiency of energy transfer as a function of the energy gaps Elhi — Erc and Elh 2 — Elhi, for the cases (a) NNN, 
(b) NNT, and (c) TTT, all in the S parameter set. The efficiency is sensitive to the site energies, with the maximum occurring 
when there is a clear energy funnel, Elh 2 > Elhi > Erc, provided that the differences in energy are not so large as to inhibit 
EET through a diminished spectral overlap. The original energy gaps are indicated by the black dots, and the features A-D are 
discussed in the text. The corresponding plots for the R parameter set are Supplementary Eig. 2. 


between LHl and the RC. 

Of course, the spectral overlap does not offer the com¬ 
plete explanation. Although the smallest energy gap 
between LHl and RC occurs when both complexes are 
trimmed, the highest LHl ^ RC rate is when both are 
delocalized. This is due both to the brightness of the 
exciton states in both complexes (supertransfer) and the 
larger number of pathways toward the RC when there 
are four BChls as opposed to two. The importance of 
the number of pathways (i.e., entropic driving) is further 
illustrated in Supplementary Table 3. 


Energy optimization 

The preceding discussion suggests that excitonic cou¬ 
plings affect EET rates more through energy-level shifts 
than through supertransfer. If so, the poor performance 
in cases such as NNT and TNT should be correctable by 
adjusting the site energies. 

We tested this hypothesis by repeating the simulations 
with site energies as free parameters. BChl site energies 
vary widely in photosynthetic complexes: single substitu¬ 
tions on the protein backbone can significantly affect the 
energy, especially if the residue axially ligates the magne¬ 
sium or if is charged and thus modifies the electrostatic 
environment m- It is therefore plausible that natural 
selection could have modified the site energies (within 
limits, of course) if it increased fitness. 

As Fig. shows, choosing the optimal value for the site 
energies of LHl and LH2 (relative to the average special- 
pair energy F^rc) dramatically enhances the efficiency in 
all cases. In the S parameters, the optimized efficiency is 
always above 93%. The R parameters do not perform as 


well due to weaker long-range couplings, but the optimal 
efficiency is nevertheless always above 73%, i.e., the same 
or greater than the natural-case efficiency of 73%. For 
details, see Supplementary Table 2. 

The behaviour of the efficiency with changes in the 
site energies F^lhi and F^lh 2 = ^(F^a + F^p) is shown in 
Fig. a Most features of the plots can be understood 
as consequences of energy funnelling. For example, the 
efficiency is negligible unless LHl is higher in energy than 
the RC or only slightly (a few a) lower, but not so high 
that their spectra no longer overlap. 

When LH2 is also considered, there are four regions of 
interest, labeled A-D in Fig. depending on the offset 
F^lh 2 — F^lhi- The peak B occurs when there is a clear 
energy funnel toward LHl, whereas in the adjacent valley 
C, LH2 is slightly lower than LHl, so that outward transfer 
LHl ^ LH2 dominates. This problem is compounded by 
the number of available pathways (i.e., entropy): because 
there are 10 LH2s in the model, an exciton on LHl can 
move to an LH2 in 10 different ways (at the same rate), 
while an exciton on LH2 has only one pathway to LHl. 
Therefore, the excitons accumulate in the LH2s, making 
the transfer to the RC unlikely. 

When the offset F^lh 2 — F^lhi is large, whether positive 
(A) or negative (D), the small spectral overlap between 
LHl and LH2 gives a small efficiency. In those cases, 
only excitons starting out in LHl or the RC contribute to 
the efficiency, while those from LH2—by far the largest 
proportion—cannot leave and are wasted. The efficiency 
in region A is higher than in region D (and increases 
further in the unphysical limit F^lh 2 ^ oc) because there 
are fewer photons at higher energies, meaning that fewer 
excitons start out at LH2 and are wasted. 

The plots for the three cases in Fig. (and for all 






















the cases that are not shown) are qualitatively the same, 
differing only in the positions of the peaks and their 
widths. The differences in peak positions are caused by 
the excitonic splittings, while the different widths occur 
because the couplings allow both lower- and higher-energy 
excitonic states to act as EET donors and acceptors. 

Role and evolution of coherence 

Delocalization in LH2 prompted speculation that it 
plays a functional role, having been selected by evolution 
because it enhanced light-harvesting efficiency. However, 
this is not so, since trimming the LH2 always increases 
the efficiency. Why then is there delocalization? The 
main evolutionary pressure is probably not on efficiency 
but on the total number of excitons processed by each 
RC. Although removing half the BChls might increase the 
efficiency marginally, it would also halve the absorption 
cross-section of LH2, decreasing the total exciton flux 
into the RC. In other words, delocalization in LH2 is a 
spandrel m, a byproduct of packing BChls densely to 
maximize absorption, just as the red color of vertebrate 
blood is not adaptive, but a byproduct of the oxygen¬ 
carrying ability of haemoglobin, which happens to be 
red. Indeed, a clade of cryptophyte algae underwent a 
mutation that reduced excitonic delocalization in their 
antennas [66] with no apparent decrease in fitness. 

By contrast, the N geometry in the RC always outper¬ 
forms T because the Davydov splitting brings its energy 
closer to that of LHl. However, it would put the cart 
before the horse to conclude that RC energy splittings 
are an adaptation to create better energy alignment with 
LHl. Rather, a strongly coupled special pair is a feature 
conserved across photosynthetic organisms [T], meaning 
that it arose before the advent of purple bacteria. 

Instead, we should ask whether purple bacteria built 
the antennas around the RC so as to create an effective 
energy funnel. The answer is yes, as confirmed by the 
high efficiency in the natural case. Even so, the alignment 
is not optimal, and Eig. [^indicates a substantial further 
increase in efficiency is possible in principle, depending on 
how much further the site energies could realistically be 
altered. Nevertheless, the continued survival of the species 
indicates that the efficiency is probably good enough 
that an additional increase would not confer a decisive 
evolutionary advantage. 

In particular, having the RC higher in energy than 
LHl creates a rate-limiting uphill step at the end of the 
EET chain. It has been proposed that the uphill transfer 
is also adaptive, preventing too many excitons arriving 
at a ‘closed’ RC (one having recently undergone charge 
separation), thus averting excessive energy dissipation 
and thermal damage [20l [211 ESj- However, we doubt 
that the exciton flux would be large enough to damage 
the RC, considering that the light-harvesting processes 


are completed orders of magnitude faster than the mean 
time between the arrivals of two photons. Therefore, it 
appears that the final uphill step is, like delocalization, 
an evolutionary relic in a light-harvesting system that 
functions well enough, but not optimally. 

CONCLUSIONS 

EET efficiency in purple bacteria can be largely under¬ 
stood in terms of the energy landscape, with supertransfer 
playing a minor role. Because ERET rates are exponen¬ 
tially sensitive to energy offsets, small changes in site 
energies—well within the range found in nature—can of¬ 
ten improve the rates more than the maximal effect of su¬ 
pertransfer. Despite this sensitivity, the purple-bacterial 
energy funnel is robust, and would yield high efficiencies 
across a considerable range of site energies. The domi¬ 
nance of energy funnelling is seen using two very different 
parameter sets, confirming that it is a sturdy conclusion 
insensitive to the details of the model. 
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Supplementary Figure 1. Energy levels of the natural (a) and trimmed (b) complexes in the R parameter set. An important 
difference with respect to the AS-parameter levels shown in Fig. 2 is that LHl ^ RC transfer is significantly uphill in the trimmed 
geometry, explaining the low efficiency of the TTT configuration in the R parameter set. 
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Supplementary Table 1. Efficiencies ( 77 ) and inter-complex energy transfer times when the site energies are the original energies 
listed in Table 1. Our results for the natural geometry NNN in the S parameter set (first line above) are consistent with those of 
Schulten’s group. In particular, our departures from Schulten’s approach—treating fewer light-harvesting complexes and using 
the simplified overlap integral in Eq. 4—are justified because they do not substantially affect the results. Our efficiency and our 
transfer times for LHl ^ RC and RC ^ LHl are similar to Schulten’s values of, respectively, 78-91% (depending on vesicle 
structure), 20ps, and 1.4ps [2T]. Our transfer times involving LH2 are shorter because we assumed a tighter packing based 
on AEM measurements (compare with Schulten’s 10ps for both LH 2 ^LH 2 and LH 2 ^LH 1 [ 2 j). As noted in the text, this 
difference largely disappears if the distances are the same. The S parameters tend to agree with experiment somewhat better 
than the R parameters, although not completely. Eor example, S parameters get the LHl ^ RC transfer time approximately 
right (estimated at 20 ps at room temperature [57l [67] L but the reverse rate (7-9 ps [ 68 ]) is much better captured by the R 
parameters. Measured transfer times involving LH2 (5ps for LH2^LH2 [69] and 3.3 ps for LH2^LH1 [70]) also agree with 
those predicted from the S parameters at the separations we used. 
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Supplementary Table 2. Efficiencies ( 77 ) and energy transfer times when the site energies are optimised for each geometry. 
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Supplementary Figure 2 . As Fig. 4, but for the R parameter set. 
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Supplementary Table 3. Effect of the number of pathways on the efficiency. As noted in the text, having more pathways going 
away from the RC than towards it implies that entropy drives charges away from the RC. To illustrate this further, we increased 
the number of BChls in the trimmed RC by replacing the ones that were discarded, indicated as T^. In other words, we used 
ordinary FRET on the natural RC geometry, even if that is not the correct description of EFT. It can be seen that doubling the 
number of BChls in each RC roughly doubles the transfer rate LHl ^ RC with respect to the values in Supplementary Table 
This results in an appreciable increase in the efficiency, averaging 75% for the S parameters and 65% for the R parameters 
(relative to Supplementary Table [^. 
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Geometry 
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Supplementary Table 4. EET efficiencies when the LH2 and/or LHl aggregates are trimmed more aggressively. denotes 
aggregates in which only every third BChl is kept and the results are shown using both the original site energies and the optimal 
ones. The efficiencies are similar to those obtained with the every-second trimming described in the text, indicating that our 
conclusions are not sensitive to the extent of trimming. 



